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ABSTRACT 


Finite-amplitude standing waves in air at ambient temper- 
atures contained within a detunable, rigid-walled rectangular 
cavity were experimentally investigated. The pressure wave- 
forms resulting from excitations in the vicinity of the reson- 
ances of the 100 and 010 modes were analyzed for harmonic 
content at three different cavity configurations. Similarly, 
waveforms of the 110 and 210 modes were analyzed for four 
different cavity configurations. These results were compared 
to the predictions of a theory of Coppens and Sanders wherein 
the resonance freguencies and quality factors of all modes were 
experimentally determined from infinitesimal-amplitude measure~ 
ments. It was found that in the case of nearly degenerate 
modes, coupling within the cavity allowed for the excitation of 
modes belonging to families other than that containing the 
driven mode. This feature is absent from the theory in its 


present form. 





Cabo gen. CONTENTS 


ies OU EUONS Rive «6 8 6 ss Sl 
adele, Pee nGRGUND AND THEORY . . « « G0 
ere EXPERIMENTAL CONSIDERATIONS ... 
i ETE SURES) 1 a a eS 
eo MECROPHONE, OUTPUT AND STRENGTH 
Pee ee OUEBNCY PARAMETER oa). «ait 
Pree oO ATION 2. 2 se 8 el el 
By. Die PAK ING (PROCEDURES . « « « « » 


Pree @imeehOCEDURES .«. . « « » « 


® ® e ® ° 


PARAMETER 


Pee ota AMPLITUDE ANALYSIS ... .« 


C. FINITE AMPLITUDE ANALYSIS .. 
We PiSWBCwMP ECs «os «& © « 6 6 3 
DEMO ONCLUSIONS . 5 . . . se ae ew 
PURO GRAPHY). . 9. 6 6 6 6 6 et ek 
Mian DESTREBUTION LIST ...... « 


FORM DbD 1473 ° ° ° e .8 ° ° e e ° » e ° r 


\O 


fe: 


14 


14 


Ais, 


L5 


a6 


18 


20 


apt 


59 


60 


61 





Pee PRODUCTION 


The purpose of this research is to investigate finite- 
amplitude standing wavesS in air at ambient temperatures in im- 
perfect rigid walled cavities for comparison with the theory 
of Coppens and Sanders [1]. This investigation was a contin- 
uation of the experimental works of Lane [2], with modifications 


mMemctake into account the non-ideal nature of the cavity. 





II. BACKGROUND AND THEORY 


A plane elastic wave traveling in a non-disSipative medium 
will change wave form as predicted by hydrodynamic equations. 
Tf the problem is extended to absorptive media, only waves of 
relatively high amplitude will change waveforms. Thus, a 
finite~amplitude wave is a solution of the nonlinear wave 
equation that describes the real life environment that con- 
fronts large amplitude acoustic waves. This finite-amplitude 
wave requires second order terms, or higher, to describe the 
waveform. 

The study of finite-amplitude effects in standing waves is 
the logical extension of earlier studies dealing with traveling 
waves [3, 4, 5, 6, and 7]. The theory of Coppens and Sanders 
[1] deals with finite~-amplitude standing waves in rigid walled 
Cavities. It is an extension of the Kechk-Beyer [8] pertur- 
bation approach which makes use of perturbation methods and 
Fourier series representations of the waveform. Initially, the 
theory included wall losses as predicted by Rayleigh~-Kirchoff 
theory [9]. This approach was investigated by Ruff [10] to 
the tenth harmonic, at which the mathematics became divergent. 
Coppens [ll] expanded the theory and used a Fourier decompo- 
Sition technique that was experimentally investigated in tubes 
by Beech [12] and Winn [13]. The tubes investigated did not 
behave in a manner predicted by Rayleigh-Kirchoff loss mechan- 
isms. The theory was revised to include empirically determined 


losses which was subsequently investigated by Lane [2]. Lane 





investigated standing waves in a rigid rectangular cavity and 
found that the theoretical model successfully predicted the 
major features of the harmonic content for finite-amplitude 
Standing waves. 

The theory of Coppens and Sanders developes a one-dimen- 
Sional, non-linear wave equation with dissipative terms which 
correspond to those encountered in the description of plane 
Standing waves in a rigid-walled cavity. The wave equation is 


then Fourier decomposed to take the form 


: — = =, Ee + Dd) “a = _ ae)” (1) 
n=l 0x Cae a ax 
where 
i = 6 (ay hn ay, ~ 2%) 
2 : om oX et AC 2 Ace 


and contains a dispersive term and a dissipative term, both 


based on the Rayleigh-Kirchoff theory of wall lesses [9], and 


64 = 1/04 

a | ieee 
a = particle displacement due to n habnonsc 
— x en 


pe= (1 + Ca Sy 


os = specific heat of medium at constant pressure 
Cc, = specific heat of medium at constant volume 
ee frequency (angular) of fundamental resonance 
QO = cee 
1 Way 7 OF 
Wy 7 YZ is the difference of the half power frequencics above 


and belcw the fundamental. 








Pamerion 1 can be reformulated [2] to introduce the actual 
Seeeity resonances and generalized to include an arbitrary, 


empirical absorption coefficient for each mode. The result is 
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where 
: ela) 
ar. = resonance frequency of the n OVeCE cone 
ee os at such that pressure antinodes exist at x = 0 and at 
> ian Wy 
Aw = nw-w 
n me 
w = frequency at which system is being driven (near Wr? 
“ny 
oo = — 
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a and Ww. are the half power frequencies above and below W ar 


The acoustic pressure is the 


co 


eee M Y R cos K_ (L-x)cos n(wt + @_) 
ox 2 n n n 
‘ee n=1 
o°O 
p = acoustic pressure 
BS Gensity of medium 
on = thermodynamic speed of sound in medium 
M = mach number = p/p ma 
EUG oni e" 
P, = peak pressure of fundamental component of wave 
Ry = FORUraecrE COCET LCIent OF aes harmonic component, normel- 
1zed such that Ry = ] 
and 
_ nae ; | : 7 
B. = phase angle of n harmonic component, where oy = 0. 





This equation is then manipulated into a general, time-inde- 
pendent form which, when evaluated at the rigid boundaries of 
the cavity, can be further simplified to produce the following 


form amenable to computer calculation: 
n=] 


HR COs (f-8) = MbQ, wae! 2 a Ra Ne eoke (B. + eee 
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fp sles 
me a Similar form with "sin" in place of "cos" 
where 
2 Aw z O 
H = BN A Bele 
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ang n 
and 
2 tin O73 
tan Ce = es as 





04/0, 
The above equations show the importance of knowing the oe and 


Wor Simthic SyVStemM LOr 1npur, tO the computer program. For 


convenience, a quantity, EA! is defined to indicate the position 
of ae relative to the classical harmonic frequencies, Wy» 
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The values of oe and Ey are directly determined from the infin- 
itesimal-amplitude analysis of the cavity and are the input 


quantities for the computer calculations. 





Pit. EXPERIMENTAL CONSIDERATIONS 


A. APPARATUS 

The rectangular cavity 1s shown in Figure 1. It was con- 
Structed with 0.985 inch thick aluminum plates, milled and 
Screwed together to provide tight, right angle joints. The 
joints were sealed with thin layers of high-vacuum silicon 
grease which was applied to all ajoining surfaces prior to 
assembly. The interior dimensions were 6.97cm by 20.96 cm by 
25.l2cm. The theoretically predicted modes assuming an ideal 
rectangular cavity are shown in Table 1. The accesses were 
provided to the cavity for the piston and microphone. These 
were positioned at opposite corners of the cavity. 

A block diagram of the experimental apparatus is shown in 
Figure 2. Acoustic waves were generated within the cavity by 
a plane-~faced piston. The piston gained access to the interior 
eeetne Cavity by way of the larga orifice. It was fitted with 
ene O-ring which was lubricated and further sealed with high 
vacuum silicon grease to provide an airtight fit. An Endeveco 
Model 2215 accclerometer was mounted on the inner face of the 
piston which enabled the piston motion to be monitored on a 
Hewlett Packard Model 400D voltmeter, a Fairchild Model 766H 
Gual trace oscilloscope, and a Hewlett Packard 302A wave 
analyzer. The piston was attached to a MB Electronics Model 
2120MB power amplifiers operating in parallel. The frequency 
and power level of the piston were controlled by a General 


Radio 1161-A coherent decade frequency synthesizer. The 
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FIGURE 1. Rigid Walled Cavity 
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fundamental and harmonic components of the pressure wave were 
monitored by a 1/4 inch diameter Bruel and Kjar Type 4136 con- 
denser microphone. It gained access to the cavity via the 
small orifice and its active element was flush with the inte- 
rior of the cavity. It was fitted with an O-ring and sealed 
with high vacuum silicon grease. The output of the microphone 
system was monitored by a Hewlett Packard Model 400D vacuum 
tube voltmeter, the Fairchild Dual Trace Oscilloscope, and four 


Hewlett Packard Model 302A wave analyzers. 


Pee Si LCROPHONE OUTPUT AND STRENGTH PARAMETER 
The investigation of the pressure waveform in the cavity 


required the calculation of the strength parameter from the 


observable quantities. The strength parameter is defined as 
ea 
M = re} 
a @ 


The microphone sensitivity was determined by piston phone cali- 





bration 
-3 2 
Sia =i 4Go82) O0s) 0 VOLE thea o> 
where 
.. £Y 
2 | en ee 
rms 


and V is the rms voltage reading on the wave analyzer. The 
strength parameter can now be reformulated with observable or 


calculable quantities as 


Spey 2 VbQ,/(S_pe*) (2) 
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fee REOUENCY PARAMETER 
The frequency parameter indicates the position of the 
driving frequency relative to the fundamental resonance fre- 


quency of the system. The frequency parameter is defined as 


FP = 2(f-nf,) Q,/nt, 


Pee CAVITY ISOLATION 

mo Manimize the possible effects that could occur due to 
Emeevyibration coupling of the cavity and exciter or cavity and 
table, the cavity and piston exciter was mechanically isolated 
with alternate slabs of wood and rubber. This restricted the 
Mmeeraction of the cavity with the piston to the piston orifice 


Say « 
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Vw lst IING PROCEDURES 


Pee RERUN PROCEDURES 

The modes of primary interest in this experiment were those 
that are in close proximity to one another, such as the 050 and 
pegeer the 330 and 420 modes. The relative position of a pair 
of modes could be adjusted by changing the position of the 
piston relative to the interior of the cavity. A quick fre- 
quency sweep at infinitesimal amplitudes revealed the relation- 
Ships between the modes of interest. This procedure was 
repeated until the desired "tuned" condition was obtained. 

> On two occasions significant changes in the frequencies of 
the modes occurred Benne Geen ee PTs was assumed to be caused 
by an unintentional movement of the piston with respect to the 
cavity. This was concluded as on both occasions the cavity 
Medebeen repositioned on the morning of the run. To correct 
this problem, the cavity was run at infinitesimal amplitudes 
for an eight hour period and then allowed to sit overnight 
prior to a data SoE Gea ALTerwanroLeating this “action, no 
further "“settling-in" problems were encountered. 

5» As with Beech [12] and Lane [2] significant drifts in 
resonant frequency were noted during each session. This is 
most likely due to temperature variations. Prior to each data 
run, the cavity was activated at infinitesimal amplitudes for 
at least four hours. This prewarmed the cavity and was found 
to significantly reduce the rate of change of the fundamental 


frequency. However, this did not preclude the necessity of 


ibs: 





taking numerous measurements of the fundamental during each 
run to establish a calibration curve. It was assumed that all 


harmonics drifted proportionally with the fundamental. 


Peeeetne LINITESIMAL AMPLITUDE ANALYSIS 

An initial graphical representation of the linear response 
of the cavity was desirable as a guide in later more detailed 
measurements. This was achieved by setting the drive level at 
eeewe 0.01 volts as indicated on the accelerometer VTVM. A 
plot was then made of the relative amplitudes (as measured on 
fhe macrophone VTVM) versus frequency. This plot gave a visual 
display of the or and En of each mode (see Figure 3). 

If this plot looked sufficiently interesting, detailed 
measurements of the or and ch for each mode were performed. 
The microphone output voltage for cach mode was maximized by 
varying the frequency. The frequency was then shifted until 
the voltage measured on the wave analyzer (in AFC) was -3dB 
below the maximum value. The difference between the frequencies 
Seeene —3d6B points allowed the calculation of Q. It was assumed 
that the center frequency between the -3dB points was a good 
estimate of the resonance frequency. In some cases where the 
resonant peak was not symmetric, measurements were taken at the 
=3€0B, -2dB and -1dB points. A linear plot of center frequency 
versus dB-down was extrapolated to 0dB to ascertain the reson- 
ance frequency. In a similar fashion a linear plot of 09 versus 
adB-down was constructed to give an estimate of Q. The equations 


msea £Or this extrapolation were 
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Lekaing numerous measurements of the fundamental during each 
run to establish a calibration curve. It was assumed that all 


harmonics drifted proportionally with the fundamental. 


Poet NE INITESIMAL AMPLITUDE ANALYSIS 

An initial graphical representation of the linear response 
of the cavity was desirable as a guide in later more detailed 
measurements. This was achieved by setting the drive level at 
about 0.01 volts as indicated on the accelerometer VTVM. A 
plot was then made of the relative amplitudes (as measured on 
the microphone VIVM) versus frequency. This plot gave a visual 
display of the = and En of each mode (see Figure 3). 

If this plot looked sufficiently interesting, detailed 
measurements of the om and Ce for each mode were performed. 
The microphone output voltage for each mode was maximized by 
varying the frequency. The frequency was then shifted until 
the voltage measured on the wave analyzer (in AFC) was —-3dB 
below the maximum value. The difference between the frequencics 
of the -3dB points allowed the calculation of Q. It was assumed 
that the center frecuency between the -3dB points was a good 
estimate of the resonance frequency. In some cases where the 
resonant peak was not symmetric, measurements were taken at the 
sob, —-2daB and -1dB points. A linear plot of center frequency 
versus dB-down was extrapolated to OdB to ascertain the reson- 
ence frequency. In a similar fashion a linear plot of 0 versus 
dB-down was constructed to give an estimate of 9. The equations 


used for this extrapolation were 
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3dB Oe = 
i - WwW 
ir n- 
Yar 
2aB Q = 0.757 ——-——-—— 
n W - WwW 
nt NG oe 
Ynr 
dB O = 07520 a Vo ee 
nN Ww = 
Olam Gls 
where 
Ww ~ 110° 0 anh ana ow © ahglelanhe a are the frequencies 
n+ n- n+ - 


corresponding to the 2dB and 1daB points above and below ie 
The En were calculated from the resonance frequencies of the 


modes. 


See ef aNITE-AMPLITUDE ANALYSIS - 

After selection of a family of modes to be investigated, 
the Q of the gravest member was measured. From this the peak 
pressure of the fundamental was calculated from Equation 2 for 
the desired strength parameter. The choice of strength para- 
meter was dependent upon the extent of the frequency sweep 
desired. As the frequency was moved off resonance, it was 
necessary to increase the drive level to keep the level of the 
Brementa 1 constant. There is a limit to how far above and 
below resonance the frequency could be changed. At high 
driving levels the accelerometer output would distort The 
onset of distrotion was abrupt and easily observed on the 
oscilloscope and was also accompanied by an easily detectable 
Seo tOLelon Of cound Coming trom the cayity. This distortion 
always occurred at drive levels below the maximum permitted by 
the power amplifier. The higher the strength parameter, the 


Narrower would be the distortion limited band inspected. 
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A bank of four wave analyzers was used allowing measurement 
Mmmmpene 1st, 2nd, 3rd and 4th or lst, 5th, 6th and 7th harmonics 
Simultaneously. This promoted consistency between amplitude 
measurements. The peak pressure was maintained constant and 
the harmonic content was recorded for a range of frequency 


about the resonance frequency of the gravest member. 
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V. RESULTS 


Figures 4 and 5 are the graphical representations of the 
cavity's linear response when excited at infinitesimal ampli- 
tudes in the Dareaneey of the theoretically degenerate modes 
feveancd 050. Cavity configurations A, B, and C correspond to 
oeerecrent piston positions within the orifice. Configuration 
A has the piston flush with the interior cavity and B has the 
Mecton withdrawn into the orifice, and in C the piston is with- 
drawn even further into the orifice. In Figure 4 each Con- 
figuration has been aligned with six times the resonance fre- 
quency of the 100 mode. Figure 5 has been aligned with five 
times the resonance frequency of the 010 mode. Observe that 
both amplitude and relative position of each peak changes sub- 
meantually with different piston positions. Significantly no 
clear distinction can be made by simple observation as to 
which family of modes each peak belongs. 

Figures 6 and 7 represent cavity response in the vacinity 
OL three times the resonance frequency of the 110 mode or two 
times the resonance frequency of the 210 mode. Each configur- 
ation has been aligned as explained in Figures 4 and 5. fn 
Ponrguratien 0, the piston is nearly flush with the interior 
of the cavity and Gonfiguration G corresponds to the furthest 
displacement of the piston within the orifice. Observe that 
with Configuration F, both peaks have merged and in Configuration 
Ga third peak has appeared. This peak is thought to be the 


040 mode, but has not been positively identified. 
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Tables II through XV present the data for the linear, 
infinitesimal-amplitude region. The information on the empiri- 
cal losses is contained in the oe and EWS Each table corre- 
Sponds to a particular family of modes. When an extrapolation 
to OdGB-has been made, 3dB, 2dB and 1dB data has been included. 
For Configurations E and F mode overlap was so complete that 
it waS necessary to estimate Q from average values of Q for the 
Same mode for other adjacent configurations. Estimation of E 
was carried out by inspection of Figures 6 and 7. 

The values from the tables were used in the computer pro- 
mam to predict the harmonic distortion curves in Figures 9 and 
11 through 23. The theoretical curves in these figures are 
plotted along with experimentally measured values of harmonic 
content. The theoretical predicted results are indicated as 
solid curves and measured values were plotted as solid symbols. 
Data were taken and theoretical predictions were made for 
strength parameters of 0.50 for driving the 100 mode, 010 mode 
and the 110 mode. A strength parameter of 0.25 was used for 
the 210 mode. 

An initial problem was to ascertain which peak belonged to 
a particular family when two or more peaks occurred in close 
proximity as seen in Figures 4, 5, 6, and 7. This was solved 
by using the 90's and E's of each peak to generate a theoretical 
prediction. The theoretical prediction was then compared to 
the experimentally measured values. If good agreement between 
predicted and measured values was found, then the peak was con- 
Sirocred 1dentitiaed. bigures S$, 9, 10, and 11 illustrate this 


technique with Configuration C. 
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The experimental data in Figures 8 and 9 are identical and 
were obtained by driving the cavity at frequencies near the 
resonance of the 100 mode. The theoretical results are cal- 
culated from the same 07 and EA exigepe LomeO(6) and EF (6). 
Peaure S contains the 0(6) and E(6) calculated from the left 
meew oO: Configuration C on Figure 4. Figure 9 contains the 
0(6) and E(6) calculated from the right peak. 

The good agreement between measured and predicted values 
exhibited by Figure 9 is clear evidence that the right-hand peak 
in Figure 4 represents the 600 mode. Additional verification 
is shown in Figures 10 and 11 which contain experimental data 
measured by driving the 010 mode and the same On and E, except 
memes) and £(5). Figure 10 contains the 0(5) and E(5) £rom 
left peak of Configuration C on Figure 4 and Figure 11 contains 
the 0(5) and E(5) from the right peak. The good agreement 
shown in Figure 10 identifies the left peak in Figure 4 as be- 
Tenging to the 050 mode. 

In a Similar manner, all peaks with the exception of Con- 
figuration G were identified and recorded on Figures 4, 5, 6, 
and 7. The E, in the case of Configuration G were so large 
that positive identification was not possible. 

Fagures 9, 10, Zee eee and) > are comparisons (Of mea 
sured data and theoretical predictions of the 100 modes and 010 
MOVeceroOn —CONUIGMNiaclOns Ayes, and C. “Figures 15°through 21 
are comparisons of the 110 modes and the 210 modes for Config- 
Uedetoncu rie he saneuG we eetmauhes 9, 10, 12, 14,15, 275. and 
20 exhibit good agreement between predicted and measured values. 


Pages lomo Geil eee and 23 (exiibit Significant differences 
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between measured data and theoretical predictions in the 
vacinity of the peak of the non-family member. This indicates 
that some coupling has occurred within the cavity which has 
populated the non-family member. Figure 19 is considered to 
be in agreement with theory, but differs due to the estimated 
value of E. Figure 23 indicates a coupling with the 210 modes 
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VI. CONCLUSIONS 


The theoretical model can be Mecdetomidentity the modes of 
Pmnon-ideal, rigid-walled cavity provided the E. Aves Stet i 
ciently small (see Figures ppm lO, and 12). 

The theoretical model in its present fori famic tOraccoun. 
for the excitation of modes other than those belonging to the 
family of the driven mode. This excitation was observed to 
occur only in the case of nearly degenerate modes. Teas 
believed to be caused by some linear coupling mechanisin within 


Pine cavity. 
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waveforms of the 110 and 210 modes were analyzed for four 
different cavity configurations. These results were compared 
BomenempLreadictilons or a theory of Coppens and Sanders 
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all modes were experimentally determined from infinitesimal- 
amplitude measurements. It was found that in the case of 
nearly degenerate modes, coupling within the cavity 

allowed for the excitation of modes belonging to families 
other than that containing the driven mode. This feature 

is absent from the theory in its present form. 
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